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Abstract. We present the results of a submillimeter dust continuum study of a molecular ridge in IC 5146 carried 
out at 850 /tm and 450 /im with SCUBA on the James Clerk Maxell Telescope (JCMT). The mapped region is 
~ 14' x 2.'5 in size (~ 2pcx0.3pc) and consists of at least four dense cores which are likely to be prestellar in 
nature. To study the dust properties of the ridge and its embedded cores, we combined the dust emission data 
with dust extinction data obtained by Lada et al. (1999) from the NIR colors of background giant stars. The ridge 
shows dust extinctions above ~ 10 mag, rising up to 35 mag in the cores. 

A map of dust temperatures, constructed from the continuum flux ratios, shows significant temperature gradients: 
we find temperatures of up to ~ 20 K in the outskirts and between the cores, and down to ~10K in the cores 
themselves. The cores themselves are almost isothermal, although their average temperatures vary between 10 K 
and 18 K. We used the extinction data to derive in addition a map of the dust emissivity parametrized by 
K 1 = ksso/kv- The average value of k! agrees well with the canonical value of Mathis (1990). We find that n' 
increases by a factor of ~ 4 from ~ 1.3 10~ 5 to ~ 5 10~ 5 when the dust temperature decreases from ~ 20 K to 
~ 12 K. A Monte Carlo simulation shows that this change is significant with regard to the estimated calibration 
uncertainties. This is consistent with models of dust evolution in prestellar cores by Ossenkopf & Henning (1994) 
which predict that grain coagulation and the formation of ices on grain surfaces in the cold, dense cloud interiors 
lead to a significant increase of the 850 fj,m dust opacity. This interpretation is furthermore supported by the 
previous detection of gas-phase depletion of CO in one of the IC5146 cores (Kramer et al. 1999). Observations of 
dust fluxes at short wavelengths are however needed to verify this result. 
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1. Introduction 

Dense, cold cores embedded in molecular clouds are the 
progenitors of protostellar cores which subsequently form 
low mass stars. Many observational studies have recently 
focused on the determination of the physical properties of 
these pre-protostellar or prestellar cores. This is because 
the initial conditions just before the collapse phase deter- 
mine the properties of the protostar. There is consequently 
considerable interest in measuring their masses and ther- 
mal balance. Since cold molecular hydrogen cannot be de- 
tected in emission, trace constituents of the cloud have to 
be studied, requiring knowledge of their abundances. Two 
approaches have generally been used. The first is via the 



low- J rotational emission lines of a rare CO isotopomer. 
The major difficulty with this technique is that CO can 
freeze out onto cold dust grains in the interiors of the 
prestellar cores thus changing significantly the CO gas- 
phase abundance from its canonical value of ~ 10 -4 (e.g. 



Bcrgin et al. 2002; 


Tafalla et al. 


2002 


; Casclli et al. 1999; 


Kramer et al. 


199E 


) and also thus ct 


tanging the thermal 
)• 


balance of the cores ( 


Goldsmith 


2001 
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Since dust is by far the most important coolant and 
constitutes about 1% of the total mass, measuring the 
dust emission has proved popular as another method. The 
recent advent of large, sensitive bolometer arrays on mil- 
limeter and submillimeter telescopes now makes it possible 
to measure the weak continuum emission even from cool 
clouds showing no signs of star formation. The surface 
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regions of prestellar cores are heated by the weak inter- 
stellar radiation field (ISRF) while cosmic rays penetrate 
easily even the densest cores. The cores therefore exhibit a 
small temperature gradient between the cold core interiors 
and the warmer surfaces. This argument was quantified by 



2. Observations and data reduction 

2.1. SCUBA data 

Wc used the Submillimcter 



g- 


Leung 


1975|) and 


Zucconi ct al. 2001; 



Evans ct al. 2001). These models show the dependence 
of the gradient on, for example, the grain material and 
on the core geometry. A comparison with observed fluxes 
should therefore allow to discriminate between these fac- 
tors. Recent observations have indeed shown that many 
prestellar cores exhibit a temperat ure gradient while oth 



ers a ppear to be isothermal (e.g. Ward-Thompson et al 
2001). 



Common User Bolometer 
Array (SCUBA, [Holland ct al. l |1999| ) on the James 
Clerk Maxwell Telescope (JCMT[J) to map IC5146. On 
September, 13th and 15th, 2000, we observed a 14' x 2.5' 
region simultaneously at 850 /im and 450 /tin. The SCUBA 
array covers a hexagonal field of view of 2.5' with 91 and 
37 pixels at 450 and 850 /im respectively To create fully- 
sampled maps, the 'jiggle' mode was used, in which the 
telescope beam is moved in a 64-position pattern by the 
secondary mirror, spending 1 s at each sample point while 
simultaneously chopping at 7.8 Hz to remove atmospheric 
emission. Since the filament is elongated in right ascension 
and smaller than about 2' in declination, we chopped in 
declination using the maximum throw of ±2f5. 

Pointing, focus, and calibration checks were done on 
NGC 7027 and Uranus at regular intervals. The zenith 
opacity was derived from simultaneous skydips at 850 and 
450 /tin. It stayed almost constant at T85o = 0.18 and 
^850 = 0.24 during the first and second night respectively, 
with a repeatability being better than 3% at 850 /im. The 
half power bcamwidths were measured on Uranus to be 
16" at 850 //m and 10" at 450 /tm. The main beam is al- 
most Gaussian at 850 /im. At 450 fj.ni the main beam ef- 
ficiency is ~ 50%, with the remaining power forming an 
els in t his p a per. T he dust coagulation model of |Ossenkopj error beam around the main beam roughly ~ 30" in size. 

We reduced the data primarily using the SCUBA User 
Reduction Fac il ity (S URF; |Jenness fe Lightfoot| ( [1998D ; 
[Jcmicss et al. ( 2002| ) ) in the standard way. The data 
were corrected for individual bolometer gains (i.e. flat- 
ficldcd), for atmospheric extinction and correlated sky 
noise, and regridded and calibrated to yield the final im- 
ages in Jybeam -1 . To improve the image quality, we de- 
convolved both maps using our own CLEAN algorithm to 
deconvolve the triple beam response of the observations: 
this final step made only small changes to the maps be- 
cause there is little emission extended by more than 150 
arcsec in declination. At 450 /im, our beam model is com- 
posed of three symmetric Gaussians sim ilar to the model 
fitted by |Hogcrhcijde fc Sandclj (|200C| ). We restored the 
cleaned images back to the original resolutions. Next, we 
smoothed the images to 30" resolution using a Gaussian 
convolving kernel. And we resampled the data onto a com- 
mon grid with 15" spacing. This allows for a pixel-by-pixcl 
comparison of the 450 /im and 850 /tm maps with the Ay 
data at the same resolution and grid. 



Another difficulty is that the grain opacity at submil- 
limcter wavelengths is poorly known and may vary due to 
grain evolution in the dense, cold core interiors. At wave- 
lengths longer than about 250 /tm, the wavelength depen- 
dence of the emissivity is customarily parametrized by a 
power-law function: k\ — tto(A/Ao) _/3 . The dust proper- 
ties are thus characterized by two parameters, the emissiv- 
ity exponent (3 and the opacity kq at a given wavelength 
Ao- In general, the opacity may not even follow a power 
strict law, and so would require more than two parameters 
to describe it; given the lack of accurate multi- frequency 
submillimcter data we do not consider more complex mod- 



& Hen ling (1994, OH94 hence on) shows that the forma- 



tion of dirty ice mantles and grain coagulation increase 



the opacity at long wavelengths (cf. 


Pollack ct al. 


1994 


Stognicnko et al. 


1995; 


Hcnning et al. 


1995 


)■ 



Some studies use the spectral energy distribution in 
the FIR and submm regime to disentangle dust tc mpera- 
ture and opacity (e.g. Ward- Thompson ct al. |2002 ). These 
observations indicate a constant spectral index of (3 ~ 2. 

Here, wc try to tackle this task by combining optical 
extinction data, obtained via NIR observati ons of back- 
ground stars obtained by Lada ct "aT| ( |l999|) , with deep 
dust emission maps at 450 and 850 /tm, obtained with 
the SCUBA bolometer on the JCMT. There are two key 
strengths of this technique. First, the submillimcter flux 
ratios in themselves are an excellent temperature diag- 
nostic in cold gas, if the dust index f3 is nearly constant, 
because the Raleigh- Jeans corrections are large at these 
frequencies ihvjkB = 32 K at 450 /tin). Second, the optical 
extinctions are an accurate tracer of dust column densi- 
ties independent of the dust temperature. Combining the 
three data sets therefore allows us to determine at each 
position the dust temperature and the grain opacity, av- 
eraged along the line of sight. 

We observed an area of 2pcx0.33pc of a filamentary 
cloud in the IC 5146 complex which consists of several 
dense starless condensations. This study thus covers a 



much larger region than the previous work of Kramer et al 



2.2. Sources of errors 

The most important source of statistical error comes from 
fluctuations in the atmospheric opacity. At 450 /tm, there 
are also possible errors from changes in dish shape and 
thus in the beam pattern as the dish temperature changes, 



(1998, 199E, Paper I and II henceforth) who studied one 
core of the IC 5146 filament. 



1 The JCMT is operated by the J AC, Hawaii, on behalf of the 
UK PPARC, the Netherlands OSR, and the Canadian NRC. 
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Fig. 1. Maps of flux densities at 850 /im and 450 /im. The 450 /jm data are smoothed to the same resolution as the 
850 fim data, both maps are at 16" resolution on a 8" grid. Contours of S$5q are 20 to 140 by 20 mJy beam -1 . Contours 
of S 45 o are 100 to 900 by 100 mJy beam- 1 . The (0,0) position is a(2000) = 2ll 1 46™32 f ?5, <5(2000) = 47?33f55"0. 



and errors when deconvolving the error beam. Finally, the 
absolute flux scales at 850 and 450 /xm are derived from 
t he flux of Ur anus, and these are uncertain to 5 — 10% 
( |Sandell|^994|) . Since this error is entirely systematic, it 
cannot account for any point to point variations in cloud 
properties, only stretching monotonically the derived tem- 
perature or k' scales. 

Uncertainties in the 450 /im maps due to the error 
beam should be small, because the CLEAN deconvolu- 
tion included the error beam pattern and because we are 
only interested in data at 30" resolution, which includes 
most of the power in the error beam. In addition, most 
of the data were taken in cool nighttime conditions when 
the dish temperature is stable. 

The best estimate of the fluctuations in atmospheric 
transmission are probably the opacities at 225 GHz mea- 
sured at the CSO every ten minutes using a tipping ra- 
diometer (Archibald ct al. 20001 ). These opacities varied 
very smoothly and by small amounts, indicating very sta- 
ble atmospheric conditions during our observations. We 
applied the opacities at 850 /im measured via SCUBA sky- 
dips to the data. The scatter of these opacities, scaled to 



225 GHz ( [Archibald et al] |2002| ) relative to the fitted opac- 
ities is less than 5% during the 11 hours of the two night 
observations. At the average CSO-zenith opacity of 0.055 
and typical airmass of 1.56, this translates into a calibra- 
tion uncertainty of ~ 10% at 450 /im. Temporal variations 



in opacity are expected to be slow, of the order of 10 min- 
utes, so the error on close by parts of the maps will be 
much smaller. 

We thus estimate that the calibration error of the flux 
density ratio S450/S85O) observed simultaneously through 
the same atmosphere, is ~ 15% and ~ 10% for the ratio 
$850 1 Am since the latter is dominated by the error of the 
flux density. In fact, the resulting dust temperatures of 
between 10 and 20 K and their small scatter (Figs. 4, 5) 
also indicate that errors cannot be much larger. 



2.3. NIR data 



The NIR observa tions of IC5146 are described in Lada 
|et alj ( |1999| , |1994| ). The determined color excess E(H—K) 
is proportional to dust column density. Following con- 
vention, the dust column density is represented in terms 
of optical extinctions Ay by converting the color ex- 
cess at each position to the optical extinction using the 
normal reddening law: Ay = 15.9E(H — K) (Ricke & 
Lcbofsky 1985) with the ratio of total to selective extinc- 
tion R = A V /E(B - V) = 3.1. Although the resulting 
map accurately reflects the distribution of dust column 
density and mass, it may not accurately predict the true 
visual extinction through the cloud. This is because grain 
growth in cold clouds can alter the reddening law at wave- 
length <C 1/Um (e.g. Fig. 2 in Mathis 1990). Conversion to 



C. Kramer et al.: Dust properties of the dark cloud IC 5146 




Fig. 2. Grey scale maps of optical extinctions overlayed with contours of flux densities at 850 /im and at 450 /im. 
Here, all three data sets were smoothed to a common resolution of 30" on a 15" grid. Extinctions vary between 3 and 
45 mag. Contours of 5*850 are 10 to 130 by lOmJy beam -1 . Contours of 5450 are 80 to 900 by 80mJy beam -1 . 



near-infrared extinctions (e.g. Ajc = 1.78E(H—K)) would 
avoid this problem. In this paper, like in Papers I and II, 
we stick to Ay and therefore compare the dust absorption 
Atg5o in IC 5146 to the absorption in the NIR of general 
interstellar grains rather than to the actual value of Ky of 
grains in IC 5146. 

3. Properties of the dust 

3.1. Maps 

The IC 5146 complex of molecular clou ds lies in Cygnus 
at a distance of 460 pc ( Lada ct al. 1999|). The part we an- 
alyze d is called Northern Streamer ( Dobashi ct al. 1992 , 
1992) for its filamentary appearance which is strikingly 



visible in the map of optical extinctions of Lada et al 
(1999). We mapped at 850 and 450 /mi the inner part of 



the ridge above about 10 mag of optical extinction cover- 
ing an area of about 2pcx0.3pc. The resulting maps of 
flux densities are shown in Fig. |] at a resolution of 0.04 pc 
(16"). The filament has a clumpy appearance. About nine 
cores are discernible at both wavelengths. Most prominent 
is the bright clump at the eastern edge of the mapped re- 
gion near (5.3',— 0.8'). Many clumps appear much more 
distinct at 850 fj.ni than at 450 /im, which is probably an 
indication that the dust is cold, with the Rayleigh- Jeans 
correction suppressing the 450 /im emission. This holds for 
example for the center core near (0',0'). There are excep- 



tions. For example the clump near (4.1',— 0.7') is much 
more distinct at 450 /im than at 850 /im which probably 
indicates a much higher fraction of warm dust. 

In Fig-H we show the same data set, but smoothed 
to a common resolution of 0.067 pc (30") in comparison 
with the optical extinctions at that resolution. Currently, 
the number of detected background stars does not allow 
to create an Ay map of still higher spatial resolution. At 
least four cores show up at all three wavelengths reach- 
ing peak extinctions between 20 and 35 mag. An excep- 
tion is the core near (0.8',— 0.5') which exhibits a peak in 
extinctions of 28 mag, while it does not show maxima in 
dust emission. Almost no background stars are found near 
this position (see Fig. 1 in Kramer ( 1999| )) but only one 
strongly reddened star. This indicates high extinctions, 
but the position of the Ay maximum may be shifted. 



The mapped region shows no signs of star-forming ac- 
tivity. Two protostcllar cores associated with CO outflows 
and IRAS point sources (Dobashi et al. 1993; Lcvreault 



1983), the FU Ori source Ell-12 with a luminosity of 25 L 



and another much weaker protostar, lie just outside of 
the SCUBA maps at offsets (8. 5', -1.2') and (-5.5', 2.47'). 
The H 11 region S125 and its associated young open cluster 
IC 5146 lie more than 1° to the east. 
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Fig. 3. (a) Scatterplot of observed ratios Sszo/Av versus S^o/Ssso- Errorbars show the estimated calibration errors. 
Four dashed lines of constant kappa-ratio and varying dust temperature are also shown. These were calculated using 
Eqs.l, 2 assuming (3 = 2. Kappa-ratios range between 1 10~ 5 and 7 10 -5 in steps of 2 10 -5 . Short perpendicular markers 
denote dust temperatures which range from 6 to 30 K in steps of 4K. (b) Scatterplot of k' versus Tdust for the mapped 
positions. Dust temperatures were derived from Eq. 1 assuming (3 = 2. The k ratios were derived from Eq. 2. The 
dashed lines delineate the temperature dependence of k' for constant Ss5o/Av=(l, 4, 7) mJy beam -1 mag -1 from left 
to right, (c) Monte Carlo simulation of two data sets (see text for details). The widths of the normal distributions are 
the estimated calibration errors, (d) The resulting scatter plot of k' versus Tdust- The large cross marks the values 
derived from the mean values of Sg^o/Ay and S450/ 'S%5$. The three contours encircle the la, 2a, and 3tr areas. 



3.2. Deriving the distribution of dust properties over 
the surface of the map 

3.2.1. Basic equations 

The dust temperature Tdust is determined by the ratio of 
flux densities per beam at 850 /tm and 450 jim which are 
optical y thin up to very high column densities: 



from the Rayleigh- Jeans law are highly significant and so 
this ratio is a strong function of both the dust temperature 
and the dust emissivity index (3. To derive the dust tem- 
perature distribution over the surface of the cloud, we as- 
sume (3 = 2 as has generally been found for dust grains in 
molecular clouds (Ward- Thompson ct al. [2002 : Johnstone 
fc Ballyj|l999|; [Bandcll et al.||1999|; |Huard et al.|[l999|; |Vissei 



'450 



/8S0\ 
S 85 o V450/ 



3+flcxp(iy K/Tdustj-l 
exp(32K/T d u S t) - 1 



et al. 1998; Goldsmith ct al. 1997) and is also consistent 



(1) 



with theory A relative error of the flux density ratio of 
10% would result in an error of Td us t of 1.4 K at 12 K. 



At the low temperatures of about 10 K prevalent in Next, we used the observed optical extinctions to de- 
dark clouds and at submillimeter wavelengths, deviations rive the distribution of the ratio of absorption and extinc- 
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tion coefficients n' = Ksm/^V over the surface of the cloud 
via: 



S, 



,sr,u 



— -Bgso (Tdust )^850 



1 K 850 

1.086 k v ' 



(2) 



i?85o(7dust) is the Planck function at 850 /im and f^sso 
is the main beam solid angle. This analysis assumes (3 — 2, 
it may thus miss any variations of (3, if present. Dust flux 
densities in the FIR would be needed to analyze the fre- 
quency dependence of the dust absorptivity in more detail. 
However, we do derive variations of k' . This approach is 
supported by the dust models of e.g. OH94 which indicate 
that grain evolution in dense, cold prestellar cores leads to 
significant changes in Kgso but leaves (3 rather constant. 

3.2.2. Scatterplots 

We use the three data sets Ssso, S45Q, of the full re- 
gion at the common resolution of 30" and for each posi- 
tion on a 15" grid. Figurc||a shows a scatterplot of the 
observed ratios S^o/Ay versus S^so/S^so- Mean and rms 
values are S$5o/Ay = 4.0 ± 1.5(38%) and S^so/S'sso = 
6.7±1.6 (24%). The scatter is thus significantly larger than 
the estimated calibration errors of 10% and 15% respec- 
tively. The reason may be variations of dust temperature 
and/or of the K-ratio k', assuming that (3 stays constant. 
Using the Ay data allows to disentangle these two vari- 
ations. We plot in Fig. |]a lines of constant k' for dust 
temperatures in the range 6 to 30 K using eqs. 1 and 2. 
These lines clearly show that the observed data with its 
uncertainty is inconsistent with a constant dust temper- 
ature and kappa-ratio. High kappa-ratios are found for 
regions of low dust temperatures. The resulting scatter- 
plot of kappa-ratio versus dust temperature (Fig.|^b) also 
shows this anti-correlation. 

Care has to be taken when interpreting this plot, since 
the scatter of points is influenced by the calibration errors 
and also by the interdependency of the two equations used 
to derive the numbers for k' and Tdust- We therefore con- 
ducted a Monte Carlo simulation to create two samples 
of data shown in Fig.||c. The rms scatter of these dis- 
tributions equals the estimated calibration errors. Their 
mean values lie at the low and high end of the observed 
data but also well within its scatter (Fig. [|a) : 5450 / Ssso = 
5.5 ±15%, S S50 /Ay = 5mJy beam" 1 mag" 1 ±10%, 
and S45o/<S , 85o = 8.0 ± 15%, Ssso/Ay = 3mJy beam -1 
mag -1 ±10%. From these distributions we derived again 
a scatterplot of Td us t and k' (Fig.||d). The two selected 
distributions are significantly different at the 3er level. We 
thus find a significant drop of dust temperatures from 20 K 
to about 12 K while the K-ratio rises at the same time from 
1.3 10~ 5 to 510 -5 . 

The average of the K-ratio is 3.3 10 -5 , comparing well 
to the canonical value of 3.7 10 -5 derived from the ex- 



tinction curve of Mathit (1990) who uses the estimate of 
Hildcbrand| (|1983) for the FIR opacity: 



«a/«v = (A/250) 13 (k 2 5o/kj) (kj/kv) and 



(3) 



Ky/nj = 3.55, K250/KJ = 0.0015 
with A = 850 /im and (3 = 2. 

3.2.3. Dependence on assumed dust emissivity 

The absolute derived values of dust temperature and k' 
are dependent on the assumed value (3 = 2, but as long as 
(3 is constant, it cannot affect the derived point-to-point 
variations in properties. The same is true if the absolute 
flux scales are in error. For example, if (3 = 1.5, the set of 
curves delineating constant n' in Fig.|]a is shifted to the 
left. The resulting temperatures arc increased. They range 
between 12 and ~ 50 K and corresponding values for k' 
range between 5 10 -5 and 0.5 10 -5 . The average value lies 
at 1.7 10 -5 , significantly lower than the canonical value. 
Such high dust temperatures are not expected for the weak 
interstellar FUV field warming IC 5146 (see next chapter). 
We thus infer that a dust emissivity index of (3 = 2 is the 
more realistic value. Of course, if the assumed fluxes of 
Uranus (which set the absolute flux scales) are in error, 
a different value of beta may also be fully consistent with 
the data. 

3.2.4. Residual emission at the OFF-positions 

The dust emission at the positions of the OFF-beam 
of the SCUBA jiggle maps may be estimated fro m the 
large s cale m aps of optical extinctions obtained by Lada 
|ct al. ( 1999 ) and shown in Fig.||. To improve slightly, we 
used a recently observed, new NIR map, covering the ex- 
tended extinction of IC 5146 (unpublished observations by 
T. Huard, J.Alves, C.Lada). The average optical extinc- 
tion in the Northern area of the OFF-beam is 3.5 mag, 
in the Southern area it is 2.8 mag. For a rough estimate 
of the corresponding dust emission at 850 and 450 /im, we 
used eqs.l, 2, and assumed an average dust temperature of 
20 K and K-ratio of 1 10 -5 which is roughly appropriate for 
the inter-core regions (cf. Fig.||). The resulting flux densi- 
ties are 8mJy beam -1 at 850 /im and 61mJy beam -1 at 
450 /im corresponding to about 6% of the peak emission 
at both wavelengths. 

The resulting error may be higher in the outskirts of 
the cores, where dust emission is rather low and thus the 
relative contribution of off-source emission may be high. 
This is e.g. the case at position (4.75'/— 0.5') between cores 
#2 and #3, where flux densities drop to 6*850 = 67mJy 
beam -1 and 5450 = 517mJy beam -1 . The resulting dust 
temperature and /t-ratio at this position are 18.4 K and 
1.45 10 -5 . The average off-source contamination is ~ 12 % 
of the measured flux densities at both frequencies when us- 
ing the above rough estimate of the OFF-source flux den- 
sities. However, the ratio S^so/Ssso is largely unaffected 
and thus is the dust temperature determination. 

We did not correct for this error which leads to a slight 
but systematic underestimate of fluxes. 
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3.2.5. Maps of Td us t and n' 



3.3. Four prestellar cores in more detail 



The resulting maps of Td us t and k' are shown in Fig.^J. 
Dust temperatures are low and lie between ~ 10 and 20 K 
for the bulk of the filament. Still higher temperatures are 
reached at a few positions at the edges of the ridge but 
are probably insignificant as discussed above. The aver- 
age and rms are (Tdust) = 16.5 ± 8 K. In the center region, 
dust temperatures are low at about 10 K and there is only 
little variation. To the east of the center, the map of dust 
temperatures shows vertical stripes: temperatures are low 
at the positions of clumps #3 and #2, while the temper- 
ature rises to above ~ 20 K in the intcrclump regions and 
at the eastern edge of clump #2. The western part of the 
map is less well characterized. 



3.2.6. Discussion 

Our finding is consistent with kv and j3 being constant, 
while the absorption coefficient Kg^o is rising with drop- 
ping temperatures. 

The colder regions arc in general associated with re- 
gions of high column densities as is seen in the map of Td us t 
(Fig. ||) and in the radial profiles of the four cores discussed 
in the next chapter (Fig. ||) . Assuming that high column 
densities correlate with high volume densities, grain co- 
agulation may be of importance. In addition, gas-phase 
molecules may freeze out on dust grains at low tempera- 
tures, altering their emissivities. Sandford & Allamandola 
( |1993| ) e.g. found a sublimation temperature of CO onto 
CO-ice of 17 K from laboratory experiments under inter- 
stellar conditions. In Paper II, we indeed found evidence 
for depletion of the CO gas-phase abundance in one of the 
IC5146 cores, which we attributed to freeze-out on grain 
surfaces. 



This view is supported by the results of the dust model 
of OH94 who show that coagulation and formation of ices 
on grains in the dense and cold cloud interiors leads to 
an increase of long wavelength absorption. Depending on 
parameters like the density or ice thickness, OH94 predict 
enhancements by upto a factor of 3.4 (Table 2 in OH94) 
while we have observed an increase by about a factor of 3.8 
(Fig. |). Grain coagulation has a larger influence than for- 
mation of ices which contributes only a factor of 1.7. Dust 
opacity spectra of icy grains displayed in OH94 (Fig. 5b, c) 
show the rise in long wavelength opacity with coagulation. 
These calculations also indicate that the emissivity index 
(3 stays almost constant for ice covered grains. 

Furthermore, Fig. 5 in OH94 shows that the slopes of 
the dust opacity spectra stay rather constant in the NIR 
while the absolute opacities in the NIR rise only slightly 
when coagulation becomes important. The constancy of 
the slopes is consistent with the NIR extinctions law an- 



The positions of four embedded cores are shown in the 
map of optical extinctions (Fig.[|). Center extinctions rise 
above 20 mag and all cores show some central symmetry. 
We therefore derived radially averaged profiles of optical 
extinction, dust temperature and k' (Fig.||). Radially av- 
eraged extinctions drop by about 50% at a distance of 1' 
corresponding to 0.13pc. Mean radial values are listed in 
Table § 

The peak in Ay near (45",— 30") was excluded since it 
exhibits an east-west orientated gradient in dust temper- 
ature and k', but no central symmetry 

3.3.1. Dust temperatures 

The average temperatures of the four cores differ signifi- 
cantly: they range between 10 K and 18 K (Table Q). Mean 
temperatures of cores #1 and #3 are only 10 and 12 K, 
much lower than the globally averaged temperature of 
16.5 K. These findings correspond well with typical dust 
and gas temperatures of starless molecular cores observed 
previously. The survey of ammonia in dense cloud cores 



alyzed by Lada ct al. (1994) who found that the ratio 
of color excess E( J — H)/E(H — K) stays constant upto 
optical extinctions of 30 mag. 



by Benson & Myers (1989) shows typical gas tempera- 
tures of between 10 and 15 K. The observational derivation 
of dust temperatures is more difficult. Often, a canonical 
dust temperature of 12 K is used for quiescent dark clouds 
( Visser et al.| [200 1| ) . Dust- and gas temperatures couple 
only at high densities of more than 10 5 enr 3 ( poldsmith 
p?50l| ; [Krugel fc Walmsley| |i~984l ) . 

Cores #1 and #4 appear to be almost isothermal. 
They show an almost constant temperature profile of 10 K 
and 16 K, respectively. The other two cores show a con- 
stant inner temperature of 15 K (core #2) and 12 K (core 
#3) for R < 40", and a rise of the temperatures by 7 
and 3 K respectively at greater radii which trace the in- 
tcrclump medium when the radially averaged optical ex- 
tinctions stay constant. This rise is also prominent in the 
map of dust temperatures (Fig.|^), in particular for the 
eastern most core #2. 

For the center region of the IC5146 filament, we 
have in addition obtained a small map (0.3pcx0.3pc, 
2.2' x 2.2') of 1.2 mm flux densities which had been pre- 
sented in Paper I. Wc derived Tdust from the observed 
Siso/Ssso ratio and k' from the independently observed 
5'i20o/^y ratio (assuming again (3 = 2). The resulting 
dust temperatures vary only between 9 and 14 K. There is 
again a systematic anti-correlation between Td us t and k': 
the dust temperature varies from 9 K and k' = 8 10~ 5 to 
14 K and k' = 1 10~ 5 . The dust temperature drops slightly 
from an average of ~ 12K at 10 mag optical extinction to 
10 K at 20 mag. This gradient is shallower than the gra- 
dient derived solely from S1200/ Av in Paper I, where a 
constant, canonical dust emissivity had been assumed. 

A natural explanation for inwardly decreasing temper- 
ature gradients is that these cores are prestellar and only 
externally heated by the weak ambient ISRF. This leads to 
cool, dense dust in the core interiors surrounded by outer, 
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Aa ['] 

Fig. 4. Grey scale maps of dust temperature 7d us t and the ratio k' = Ks5o/ K v together with contours of optical 
extinctions. Contours of Ay start at 5 mag in steps of 5 mag to 45 mag. The spatial resolution of all data is 30". 
Arrows denote the center positions of four prestellar cores which are studied in more detail. Circles of 0.067 pc (30") 
radii are centered at the core positions and indicate the region used to calculate core masses. The dust temperature 
derived from the ratio of flux densities at 450 and 850 /zm is in Kelvin. The ratio k' was determined using in addition 
the optical extinction data is given in units of 10~ 5 . 

Table 1. Columns 3 and 4 give average dust temperatures and emissivities of the four cores and of the total region 
mapped. Core averages are derived from the radial profiles for R < 60". Average values were weighted with the number 
of pixels per annulus. The last line shows average values of the total mapped region covering all pixels. 
Columns 5 to 7 give core masses within a radius of 30" and masses of the total region mapped at 850 and 450 /mi. 
Mav is the mass calculated from the Ay data. M canon is calculated independently from the Ssso flux densities using 
Tdust = 12 K and Kg^o = 0.01 cm 2 g -1 . M core is the mass calculated using the derived radial profiles of Td us t and k' . 
Mass ratio are presented in Col. 8 and 9: R\ = M canon /MA v , R2 = M cme /M^- 
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12.6 ±0.4 


5.5 ±0.7 


4.97 
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17.7 
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16.5 ±8 


3.3 ±2.6 


106.5 


124.4 









less dense layers of warmer dust. The FUV field impinging 
011 IC 5146 is probab ly slightly enhanced due to the small 
IC5146 star cluster ( Wilking ct al. 1984 ) which lies at a 
distance of ~ 1° to the east. Assuming that BD±46.3474, 
a B0V star and the most massive member of the clus- 
ter, is the only source of FUV photons, that the star is 



a black body at the effective te mperature of 30 900 K cor- 
responding to its spectral type ( Panagia 1973 ), and that 
its radiation is not diluted by material along the way, we 
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Fig. 5. The left 2x2 plots show radial profiles of dust temperature (dotted lines) and optical extinction (drawn lines) 
for four cores showing peak extinctions greater than 20 mag. The right 2x2 plots show radial profiles of k' (dotted 
lines) versus optical extinction (drawn lines) for the four cores. The error bars show the rms error within each annulus, 
which includes noise and any deviations from circular symmetry in the cores. 



derive a flux of \ = 3 in units of the average interstellar 



radiatim field at a distance of 8.1 pea 



Ward-Thompson ct al. ( 2002| ) observed 18 prestellar 
cores using ISOPHOT data at 170 and 200 /mi, From the 
flux ratio, assuming (3 = 2, they also find similarly low 
dust temperatures of between ~ 10 and 20 K. About half 
of their sources show a uniform temperature distribution 
while the other half shows an inwardly decreasing temper- 
ature profile. 

Wc find low l.o. s. averaged dust temperatures of ~10 K 
in some cloud cores . This corresponds well with models of 
( 1983 ) who find that deep inside clouds, both 



Mathis et al 



silicate and graphite grains attain local temperatures be- 
tween 5 and 7K. We have found indications for increased 
opacities in the cloud interiors due to grain coagulation 
and formation of ices. This increased emissivity may lead 
to more efficient cooling of the dust grains, thus further 
lowering the dust temperatures. 



R e cent m odels by Zucconi et al.| (|200l| ) and Evans 
ct al. ( [2001 ) quantify the temperature gradient expected 



for prestellar cores illuminated only by the external av- 
erage ISRF. These models show that local temperatures 
vary between ~8K in the interior and 14 K at the surface. 
An increase of the ISRF by a factor of 2 would lead to an 
increased dust temperature by only 15-20%. 

3.3.2. Dust emissivity 

The average dust temperatures and k' values (Table 0) are 
anticorrclatcd, as had already been found when analyzing 
all the mapped positions (Fig.|J). Cores #2 and #4 show 
a low average ratio of 2 10 -5 (and at the same time rela- 
tively high temperatures) while cores #1 and #3 show a 
high ratio of 5 and 6 10 -5 respectively (and relatively low 
temperatures). The radial profiles (Fig.||) are mostly flat 
(#1, #2, #4). Only core #3 shows a large radial variation 
of the kappa-ratio between 3 and 8 10 -5 . 



The integrated flux % °f the F UV field is given here i n 

-3 



units of $ D = 2.7 10 3 erg cm 2 s 1 (Draine & Bertold 



199£ 



This unit is a factor of 1.71 stronger than the Habing field Go 
(|Habing|[i968|). 



3.3.3. 



lasses 



We calculated the masses of the prestellar cores (Fig. |], 
Table pi) by integrating over a disk with a radius of 0.067 pc 
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(30"), i.e. covering the inner region of the radial profiles 
shown in Fig.|H| We use three different methods to de- 
rive core masses. (1.) Optical extinctions were used to 
derive masses Ma v via the canonical ratio of N(H 2 )/ Ay = 
9.36 1 20 cm- 2 mag- 1 observed by |Bohlin ct al.| ( [1978D in 
the diffuse interstellar medium: AIa v = N (K2) A /j,mn 
with the area A and the mean particle mass /1 = 2.33 in 
units of the hydrogen atom mass ran- (2.) The flux den- 
sity at 850 /im was also used to independently derive the 
total mass M canon = F 85a d 2 / {K 8 5oB 8 5o{Tdust)) , where d is 
the distance, and F S50 = Ssso/^sso x A is the integrated 
flux. Here, we used canonical values for the dust temper- 
ature and k: Td us t = 12 K and k 85 o = 0.01cm 2 g _1 . The 
masses derived from dust emission are a sensitive function 
of dust temperatures. For example, assuming a tempera- 
ture of 10 K would result in an overestimate by a factor 
of 2.1 relative to the mass at 15 K temperature. Assuming 
on the other hand a temperature of 20 K would lead to an 
underestimate by a factor of 0.6 relative to the mass at 
15 K. (3.) We therefore employed a third method, using 
the radial profiles of Tdust (Fig. ||) and k' to derive masses 
M core |. 

Core masses lie between 3M Q and 7M Q , thus consti- 
tuting only about 20% of the total mass. The masses de- 
rived by the three methods agree to within 30%. This ex- 
cellent agreement is surprising, given that the masses Ma v 
and Mcanon were derived from two independent data sets 
i.e. the Ay data on the one hand and the submillimetric 
flux densities on the other hand and given the observa- 
tional errors of both methods. The agreement shows that 
the observational errors are still smaller and the underly- 
ing assumptions are valid to within that accuracy. 

When comparing Ma v with core masses M core derived 
from the profiles of dust temperature and k', we find that 
the latter are smaller by an almost constant factor of 
0.85 ± 0.01. The scatter is significantly smaller than the 
scatter of M canon / Mxv The main factor contributing to 
the improved mass estimate is probably the determined 
dust temperature. A correction of the canonical conver- 
sion factor N(H2)/j4y or, alternatively, of the extinction 
coefficient ny by only 15% would be needed to explain the 
constant deviation of A/a v from M corc . 



4. Summary 

We mapped the submillimeter dust emission at 850 and 
450 /im and the optical extinction of the IC 5146 molecular 
cloud filament. All data was smoothed to the resolution 
of the Ay data, i.e. to 0.067 pc (30"). The map covers a 
region of ~ 2pcx0.3pc. Several peaks of Ay > 20mag 
show up along the east-west orientated dust ridge which 
we interpret as dense prestellar cores. 



3 To calculate ksso from the k' profile, we assumed a constant 
extinction coefficient in the V band of Ky = 273.6 cm 2 g ~ 1 . This 
value was in turn derived using the canonical values of 
( |!990j ) (Eq.|) and a mean k 8 50 = 0.01cm 2 g" 1 . 



— The maps of dust emission correspond well with the 
map of dust extinction. Most cores show up in all three 
maps. However, the map of dust extinction shows a 
higher contrast. The subsequent analysis relies on the 
estimated calibration accuracy, i.e. mainly on the as- 
sumption that the ratio of flux densities S450 / S 8 5o is 
accurate to within ~ 15%. This accuracy is estab- 
lished by accurate correction for the atmospheric opac- 
ity, by deconvolving the error beam at 450 /xm, and by 
smoothing all three data sets to 30" resolution. The 
validity of the estimated calibration error is indicated 
by the smooth variation of dust temperatures over the 
surface of the dust ridge (Fig. 4) and the small sam- 
pling errors of Td us t within the annuli of the radial core 
profiles (Fig. 5). These errors are less than IK for all 
annuli. 

— We constructed a map of dust temperature from the 
dust cmissivitics, assuming a constant dust emissivity 
index of (3 = 2 and only one dust temperature compo- 
nent along each line of sight. A Monte Carlo simulation 
shows that the significant part of the dust tempera- 
ture distribution is restricted to temperatures between 
~ 12 K and ~ 20 K. Cores have rather low tempera- 
tures of less than ~ 15 K, while the inter clump dust 
shows temperatures of about ~ 20 K. 

— A map of the dust absorption ratio k' = k-85o/ k v was 
constructed by using in addition the map of optical 
extinctions. Its mean value corresponds well to the ra- 
tio found by [Mathis| fll990j ). The Monte Carlo simu- 
lation indicates a systematic and significant rise of k' 
with dropping dust temperature, i.e. from ~ 1.310 



at - 20 K to ~ 510 



-5 



at 



12K. 



Mathis 



This translates into a corresponding variation of kssoj 
assuming that ny stays constant. Such a variation is 
in fact predicted by the grain evolutionary model of 
OH94 when coagulation and formation of ices on grain 
surfaces becomes important in the cold, dense cloud 
interiors. 

This is also consistent with the previously observed 
reduction of the CO gas phase abundance due to freeze 
out onto grains in the IC 5146 center region (Paper II). 

— We identified four cores with high optical extinctions 
along the ridge. Radial profiles of dust temperature 
show almost isothermal profiles for the inner radius 
of 0.067 pc (30"). Two of the cores show a significant 
increase of dust temperature at larger radii which we 
attribute to the presence of externally heated dust of 
the inter clump medium. 

— Core masses Ma v derived from Ay and integrated over 
a diameter of 0.13pc (1') vary between 4 and 7M Q . 
Masses M canon derived independently from the 850 fim 
dust emissivity assuming canonical values for k 8 5o and 
Tdust agree with Ma v by better than 30%. This indi- 
cates that both completely independent methods trace 
the distribution of total dust and gas column densities 
with that accuracy. Masses M core derived from S 8 5o 
taking into account the radial profiles of Tdust and k' 
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are offset from Ma v by an almost constant factor of 
0.85. 

Dust fluxes at shorter wavelengths than ~ 450 yum are 
however needed to better disentangle variations of dust 
temperature, dust emissivity index (3, and the absolute 
opacities. This is because the spectral energy distributions 
of cold cloud cores peak around A ~ 100 — 300 fira. Since 
this wavelength range is inaccessible from the ground, it 
will be space-borne instruments like PACS and SPIRE on 
board the Herschel space observatory, which is scheduled 
to start in 2007, which will provide unprecedented sen- 
sitivity and spatial resolution to access this wavelength 



regime (e.g. Andre 2001) 
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